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Yellow mosaic disease (YMD) is one of the major diseases affecting mungbean (Vigna radiata (L.) Wilczek). 
In this study, we report the mapping of the quantitative trait locus (QTL) for mungbean yellow mosaic India 
virus (MYMIV) resistance in mungbean. An F 8 recombinant inbred line (RIL) mapping population was gen- 
erated in Thailand from a cross between NM10-12-1 (MYMIV resistance) and KPS2 (MYMIV susceptible). 
One hundred and twenty-two RILs and their parents were evaluated for MYMIV resistance in infested fields 
in India and Pakistan. A genetic linkage map was developed for the RIL population using simple sequence 
repeat (SSR) markers. Composite interval mapping identified five QTLs for MYMIV resistance: three QTLs 
for India (qYMIVl, qYMIVl and qYMIV3) and two QTLs for Pakistan (qYMIV4 and qYMIVS). qYMIVl, 
qYMIV2, qYMIVS, qYM!V4 and qYMIV5 explained 9.33%, 10.61%, 12.55%, 21.93% and 6.24% of variation 
in disease responses, respectively. qYMIVl and qYMIV4 appeared to be the same locus and were common 
to a major QTL for MYMIV resistance in India identified previously using a different resistant mungbean. 
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Introduction 

Mungbean (Vigna radiata (L.) Wilczek) is an important le- 
gume crop of South and Southeast Asia. The total planting 
area of mungbean is about 6 M ha. Mungbean is relatively 
tolerant to drought and can be harvested within 60-75 days, 
so it can be grown as a sole crop or as a component in vari- 
ous cropping systems for its dry seeds, which are rich in pro- 
teins, amino acids and minerals. Several foods are made 
from mungbean seeds. Apart from the seeds, mungbean 
sprouts are very popular in Asian cuisine. 

One of the most destructive and important diseases of 
mungbean is yellow mosaic disease (YMD). The disease 
mainly occurs in South Asian countries, viz. India, Sri 
Lanka, Pakistan and Bangladesh. The disease can reduce 
mungbean seed yield by up to 100% or even kill a plant 
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infected at an early vegetative stage. YMD is caused by 
geminivirus (genus Begomovirus, family Geminiviridae), 
which has bipartite genomes (DNA A and DNA B). Begomo- 
viruses cause mungbean yellow mosaic virus (MYMV), 
mungbean yellow mosaic India virus (MYMIV), Dolichos 
yellow mosaic virus (DoYMV) and horsegram yellow mosa- 
ic virus (HgYMV), which have spread throughout South 
Asia. MYMV and MYIMV are the main pathogens causing 
YMD in mungbean grown in India (Hussain et al. 2004, Pant 
et al. 2001), while only MYIMV is responsible for YMD in 
Pakistan (Ilyas et al. 2010). The viruses are transmitted by 
whitefly (Bemisia tabaci) and can infect mungbean at all 
growth stages. Disease symptoms after virus infection de- 
pend on the host species and susceptibility of each plant. 
Disease symptoms vary from a few small yellow specks or 
spots on a few leaves, to yellowing or chlorosis of all leaves 
of the whole plant followed by necrosis. In highly suscepti- 
ble plants, symptoms includes shortening of internodes, 
severe stunting of plants with no yield or few flowers and 
deformed pods producing small, immature and shriveled 
seeds. 
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Chemical control of the virus vector, whitefly, is ineffec- 
tive and is detrimental to the environment. Using resistant 
cultivars is the best option to manage YMD. Improved resis- 
tance to YMD is now the major goal in mungbean breeding 
programs in several mungbean production countries. Resis- 
tance to YMD in mungbean was reported to be controlled by 
a single major recessive gene with modifiers (Malik et al. 
1986, Thakur et al. 1 997), or two recessive genes (Dhole and 
Reddy 2012) or complementary recessive genes (Shukla and 
Pandya 1985). However, a major difficulty in breeding 
YMD-resistant mungbean is field screening for the disease, 
which is hampered by non-uniform development of the dis- 
ease due to fluctuation of the whitefly population in different 
locations and different seasons. Marker-assisted selection 
(MAS) is a useful tool for genetic improvement of crops. 
Selection based on molecular markers associated with the 
target trait helps in decreasing the number of phenotypic 
evaluations and thus reducing time and cost and increasing 
gain from selection. So far, there are few reports on DNA 
markers associated with YMD resistance in mungbean or 
blackgram, a related species of mungbean and all except one 
were based on gene tagging. Maiti et al. (201 1) reported that 
resistance gene analog (RGA) markers, YR4 and CYR1, 
were associated with resistance to MYMIV in blackgram 
(Vigna mungo (L.) Hepper). CYR1 cosegregated with the 
resistance in F 2 and F 3 populations from the cross between 
resistant and susceptible blackgram and was completely 
linked with the MYMIV-resistant blackgram germplasm. 
They further demonstrated that marker CYR1 was also asso- 
ciated with resistance in the mungbean germplasm. Recent- 
ly, Dhole and Reddy (2013) reported that SCAR marker 
MYMVR-583 is associated with MYMV resistance in 
mungbean. The marker is 6.8 cM from the resistance gene. 
Gupta et al. (2013) reported that CEDG180 is linked to 
MYMIV resistance in blackgram at a distance of 12.9 cM. 
In a recent study, which is the only QTL mapping study of 
YMD resistance in mungbean, Chen et al. (2013) identified 
four major QTLs on three different linkage groups for 
MYMIV resistance using AFLP and SSR markers. 

We report quantitative trait locus (QTL) mapping of 
MYMIV resistance in mungbean using SSR markers from 
mungbean and other related legumes in this paper. The ob- 
jective of this study was to locate and compare QTLs for 
MYMIV resistance in India and Pakistan. 

Materials and Methods 

Plant materials and DNA extraction 

A population of 122 recombinant inbred lines (RILs) was 
developed from a cross between Kamphaeng Saen 2 (KPS2) 
(female parent) and NM10-12-1 (male parent). KPS2 is a 
popular high-yielding mungbean variety grown in Thailand 
but is susceptible to MYMIV, while NM10-12-1 is a breed- 
ing line from Pakistan resistant to MYMIV. The RILs were 
developed at Kasetsart University, Kamphaeng Saen, 
Thailand by the single seed descent method from individual 



F 2 plants until F 8 generation. This population is the same as 
that used by Prathet et al. (2012) for mapping QTLs confer- 
ring resistance to iron-deficient chlorosis caused by calcare- 
ous soil. 

MYMIV evaluation 

YMD resistance was evaluated in India and Pakistan. In 
India, the RIL population and parents were evaluated for 
resistance under field conditions with a high population of 
whitefly and severe incidence of MYMIV at Punjab Agri- 
cultural University, Ludhiana, Punjab, India during August 
to October 2012. This location is the same site as that used 
for MYMIV evaluation, reported recently by Chen et al. 
(2013). The RILs and parents were sown in a randomized 
complete block design (RCBD) with three replications. Each 
replicate had 1 5 plants in a row of 3 m in length with 40-cm 
inter-row distance. A local susceptible check variety was 
sown after every four rows of the tested entries. Thirty days 
after planting (DAP), infection and disease severity were 
visually scored for MYMIV. In each replication, 10 to 15 
plants (depending on germination) of each entry were scored 
separately by two trained staff. The plants were scored using 
the following scale of 1 (highly resistant) to 5 (highly sus- 
ceptible) (Table 1). 

For YMD evaluation in Pakistan, the parents and the RIL 
population were evaluated for resistance at the Nuclear Insti- 
tute for Agricultural Biology, Faisalabad, Pakistan in July 
2008. They were evaluated under field conditions with a 
high incidence of whiteflies and MYMIV. Individual RIL 
and each parent were arranged in a RCBD with three repli- 
cations. Each replication had 24 to 28 plants in a row of 4 m 
in length with 40 cm inter-row distance. The parents were 
sown as checks after every two tested entries. Infection and 
disease severity were visually scored for MYMIV reactions 
at 30 DAP. The scoring system was different from that used 
in India. Experimental units (each plant of each RIL) were 
scored using scale of 1 (immune) to 5 (highly susceptible) 
(Table 1). Then, the score was used to calculate the percent- 
age disease index for each line following Akhkar et al. 
(2009). 

Genetic linkage map and analysis of new markers 

A genetic linkage map previously developed for the same 
RIL population to locate QTLs controlling resistance to cal- 
careous soil (Prathet et al. 2012) was used in this study. In 
brief, the population was analyzed with 62 polymorphic SSR 
markers screened from 1,203 SSR markers derived from 
mungbean, azuki bean {Vigna angularis), cowpea {Vigna 
unguiculatd) and common bean {Phaseolus vulgaris), to- 
gether with 8 polymorphic AFLP markers generated from 4 
primer combinations. The linkage map was constructed us- 
ing software JoinMap 3.0 (Van Ooijen and Voorrips 2001). 
The polymorphic markers were assigned to linkage groups 
(LGs) using a minimum logarithm of an odds (LOD) score 
of 4.0 and maximum recombination frequency (r) of 0.5. 
The Kosambi mapping function (Kosambi 1944) was used 
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Table 1. 


Disease scales for scoring of mungbean yellow mosaic India virus in India and Pakistan 




India 


Score 


Symptoms 


Disease reaction 


0 


No visible symptoms on leaves or very minute yellow specks (1-2) on leaves 


Highly resistant 


1 


Small yellow specks with restricted spread covering up to 5 % leaf area 


Resistant 


2 


Yellow mottling covering up to 15% leaf area 


Moderately resistant 


3 


Yellow mottling and discoloration of up to 30% leaf area 


Moderately susceptible 


4 


Pronounced yellow mottling and discoloration of leaves and pods, reduction in leaf size and stunting of 


Susceptible 




plants covering up to 75% of foliage 




5 


Severe yellow discoloration of leaves covering >75% of foliage, stunting of plants and reduction in pod size 


Highly susceptible 



Pakistan 



Score 


Symptoms 


Percentage of 
disease index" 


Disease reaction 


0 


Complete absence of symptoms 


0 


Immune 


1 


Few small yellow specks or spots on few leaves seen 


0.01-10 


Highly resistant 


2 


Bright yellow specks or spots common on leaves, easily observed and some 
coalesced 


10.01-25 


Resistant 


3 


Mostly coalesced bright yellow specks or spots common on leaves, but no or minor 
reduction in yield 


25.01^10 


Moderately resistant 


4 


Plants showing coalesced bright yellow specks or spots on all leaves, with no or 
minor stunting and set fewer normal pods 


40.01-60 


Susceptible 


5 


Yellowing or chlorosis of all leaves on whole plant followed by necrosis, shortening 
of internode, severe stunning of plants with no yield or few flowers and produce 
deformed pods 


>60.01 


Highly susceptible 



" Calculated as (sum of all disease ratings/total number of plants scored) x 20. 



to calculate the genetic distance between markers. Linkage 
groups were named following the azuki bean linkage map 
(Han et al. 2005). The analysis assigned 54 polymorphic 
markers to 12 linkage groups (LG) with a total genetic 
distance of 991.98 cM, while 16 markers were unlinked 
(Prathet et al. 2012). The shortest and the longest LGs were 
23.66 cM and 139.45 cM, respectively, with the average 
total distance per LG of 82.66 cM. The number of markers 
per LG was 2 (LG 9B) to 8 (LG 5), with an average of 4.5. 

To confirm the association between the markers, CYR1, 
YR4 and MYMIV resistance in mungbean (Maiti et al. 
2011), these markers were screened for polymorphism be- 
tween the parents. Total genomic DNA was isolated from 
young leaves of NM10-12-1 and KPS2 using the CTAB 
method described by Lodhi et al. (1994) with a slight modi- 
fication. The concentration of DNA was determined by com- 
paring against lamda DNA and diluted to 1 ng/ul for marker 
analysis. PCR was performed using a total volume of 10 ul 
containing 1 ng genomic DNA, lx Tag buffer, 2 mM MgCl 2 , 
0.2 mM dNTPs, 1 U Taq DNA polymerase (Fermentas) and 
5 pmol each of forward and reverse primers. The DNA was 
amplified in a GeneAmp PCR System 9700 (Applied 
Biosystems) under the following conditions: 94°C for 2 min 
followed by 35 cycles of 94°C for 30 s, 55°C for 30 s, 72°C 
for 1 min, with a final extension step of 72°C for 10 min. The 
PCR products were separated on 5% denaturing polyacryl- 
amide gel and visualized by silver staining. 



QTL analysis 

Markers associated with QTLs for MYMIV resistance 
were determined by single regression analysis at P = 0.001 
using i?-program 2.10.0. To locate the position of the QTLs 
on the linkage map, composite interval mapping (CIM) was 
performed with software WinQTL Cartographer 2.5 (Wang 
et al. 2010). Significant LOD thresholds value for declaring 
QTLs were computed by a 1,000 run of permutation tests at 
J P = 0.05,P = 0.01 andP = 0.001. 

Results 

MYMIV resistance in the RIL population in India and 
Pakistan 

Resistance to MYMIV in the RILs and parents was eval- 
uated under field conditions in India and Pakistan using dif- 
ferent scoring systems. In India, 1 1 lines did not germinate 
and one of the three replicates had to be discarded due to 
poor stand from soil-borne disease infection in several plots. 
KPS2 and NM10-12-1 showed highly susceptible and resis- 
tant reactions with scores of 5.00 and 1.17, respectively. The 
disease scores among the RILs varied from 0.71 to 5.00 with 
an average of 3.73 (Supplemental Table 1). The correlation 
coefficient (r) of the disease score between the two replica- 
tions was moderate, being 0.64 and significant (P< 0.0001, 
df = 108). Frequency distribution of the disease scores in the 
RILs showed skewness towards the susceptible parent 
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Fig. 1. Frequency distribution of disease scores and percentages of the 
disease index in the mungbean RIL population derived from cross 
KPS2 x NM10-12-1 evaluated for resistance to mungbean yellow mo- 
saic India virus in 2012 in Punjab, India (A) and in 2008 in Faisalabad, 
Pakistan (B). 

(Fig. 1A). Of the 111 RJLs, none were highly resistant, 5 
were resistant, 20 were moderately resistant, 9 were moder- 
ately susceptible, 42 were susceptible and 35 were highly 
susceptible. 

For the evaluation in Pakistan, KPS2 and NM10-12-1 
were highly susceptible and highly resistant with percentage 
disease indices of 80.80% and 23.64%, respectively (Sup- 
plemental Table 1). The RILs showed percentage disease 
indices between 31.89% and 100%, with an average of 
74.11%. Correlation of the percentage disease index be- 
tween replications was very high and significant (r = 0.97 to 
0.98, P< 0.0001, df=120). Frequency distribution of the 
percentage of disease index in the population demonstrated 
skewness towards KPS2 (Fig. IB). 

There was a significant correlation between disease re- 
sponses in Pakistan and India (r = 0.37, df=109, P = 
0.0001). 



QTLs conditioning MYMIV resistance 

QTL analysis by single marker analysis revealed that 
only one marker, cp02662 on LG2, was associated with 
MYMIV resistance in India. All seven markers on LG2, 
except CEDG244, were associated with MYMIV resistance 
in Pakistan. The coefficient of determination (R 2 ) of the sig- 
nificant marker in India was 8.37%. R 2 of the significant 
markers in Pakistan ranged from 9.88% for CEDG108 to 
25.10% for cp037 15. 

CIM was employed to locate QTL for MYMIV resistance 
on the linkage map. The LOD score threshold for MYMIV 
resistance in India and Pakistan was determined by a permu- 
tation test. Three QTLs were identified for the MYMIV in 
India (Fig. 2 and Table 2). The QTLs were each on LGs 2, 4 
and 9 A and were designated as qMYMIVl, qMYMIV2 and 
qMYMIV3, respectively. qMYMIVl was located between 
markers CEDG100 and cp02662, showing an LOD score of 
2.62, which was significant at P = 0.05. This QTL had an 
additive effect of 0.35 and explained 9.33% of the variation 
in disease reaction. qMYMIV2 was localized between mark- 
ers DMB-SSR008 and VR113, showing LOD score of 2.54, 
which was significant at P- 0.05. This QTL showed an ad- 
ditive effect of -0.50 and accounted for 10.67% of the vari- 
ation in the disease score. qMYMIV3 was located between 
markers CEDG166 and CEDG304 and showed significance 
at P = 0.01 with an LOD score of 3.42. It had an additive 
effect of 0.42 and accounted for 12.55% of the disease score 
variation. qMYMIVl and qMYMIV3 alleles from NM10-12- 
1 increased the resistance, while qMYMIV2 alleles from 
NM10-12-1 decreased the resistance. 

Two QTLs, qMYMIV4 and qMYMIV5, were detected by 
CIM for disease resistance in Pakistan. qMYMIV4 had an 
LOD score of 10.00 and was significant at P= 0.001. It is 
located on LG2 between markers CEDG100 and cp02662 
(Fig. 2). qMYMIV4 had an additive effect of 6.41% and 
explained 27.93% of the variation in the percentages of the 
disease index. qMYMIVS was detected on LG6 between 
markers CEDG121 and CEDG191. This QTL showed a 
LOD score of 2.55 and was significant at P = 0.05. It had an 
additive effect of 3.04% and explained 6.24% of the percent- 
age of the disease index. At both QTLs, alleles from NM10- 
12-1 increased resistance at this locus. It is important to note 
that another QTL detected on LG9A almost reached signifi- 
cance at P = 0.05 (LOD threshold = 2.35) (Supplemental 
Fig. 1). This locus was located between markers CEDG166 
and CEDG304 with an LOD score of 2.27. 

Discussion 

In this study, we located QTLs for resistance to MYMIV in 
two different countries, India and Pakistan. CIM identified 
five QTLs for resistance: three QTLs in India (qMYMIVl , 
qMYMIV2 and qMYMIVS) and two QTLs in Pakistan 
(qMYMIV4 and qMYMIVS) (Table 2 and Fig. 2). Although 
the disease evaluation was conducted in different locations 
and years and with different scoring systems, a common 
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Fig. 2. An SSR-based linkage map of the mungbean RIL population derived from the cross KPS2 x NM10-12-1 and position of the QTLs condi- 
tioning yellow mosaic India virus resistance, as detected by composite interval mapping. 



Table 2. Quantitative trait loci conditioning resistance to mungbean yellow mosaic India virus in F g mungbean RIL populations from KPS2 and 
NM10-12-1 in India and Pakistan 



Location 


QTL name 


Linkage group 


Marker interval 


Position (cM)° 


LOD score 


PVE (%) b 


Additive effect 


India 


qYMIVl 


2 


CEDG100-cp02662 


38.8 


2.62 


9.33 


0.35 




qYMIV2 


4 


DMB-SSR008-VR113 


51.8 


2.54 


10.67 


-0.50 




qYMIV3 


9A 


CEDG166-CEDG304 


50.2 


3.42 


12.55 


0.42 


Pakistan 


qYMIV4 


2 


CEDG100-cp02662 


29.7 


10.00 


27.93 


6.41 




qYMIVS 


6 


CEDG121-CEDG191 


0.1 


2.55 


6.24 


3.04 



" Position on the linkage group. 

b Percentage of phenotypic variance explained by the QTL. 



QTL was identified. qMYMIVl detected in India and 
qMYMIV4 identified in Pakistan are in the same marker in- 
terval and can be considered as the same QTL. In addition, 
qMYMIV3 detected in India showed the highest genetic ef- 
fect and was almost significant in Pakistan. Therefore, two 
out of the five QTLs detected in India and Pakistan appeared 
to be common. Moreover, qMYMIVl and qMYMIV4 detect- 
ed in our study and the major QTL reported by Chen et al. 
(2013) appeared to be the same, although only one SSR 
marker (DMB-SSR151) was common between the two ge- 
netic linkage maps. The location for evaluation for the resis- 



tance in India in our study (Punjab Agricultural University, 
Ludhiana) was the same as that reported by Chen et al. 
(2013). Based on the mungbean genetic map reported by 
Kajonphol et al. (2012), SSR marker DMB-SSR158, which 
links to the major QTL controlling MYMIV resistance, 
MYMIV r 9_25, identified by Chen et al (2013), was on LG2 
between markers CEDG100 and CEDG108. Our linkage 
groups and those of Kajonphol et al. (2012) are the same and 
CEDG100 and CEDG108 were located on LG2 of both 
maps. A QTL for the resistance in both India and Pakistan 
was located between the two markers in our study. 
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Therefore, MYMIV9_25, qMYMIVl and qMYMIV4 are 
common QTLs conferring MYMIV resistance. This locus 
was detected in different locations and years, with different 
sources of resistance and scoring systems; thus, it is highly 
robust and can be used in marker-assisted selection. DMB- 
SSR1 51, which is linked to MYMIV r 8_48.8, a QTL control- 
ling MYMIV identified by Chen et al. (2013), was mapped 
to LG5 of our map, but no QTL was detected in this genome 
region in our study. The resistance parent used in our study, 
NM10-12-1 and that of Chen et al. (2013), NM92 (a deriva- 
tive of NM36 x VC2768B), were both developed by the Nu- 
clear Institute for Agriculture and Biology (NIAB), Pakistan 
(Ali et al. 1997). NM36 is derived from a cross between ac- 
cession 6601 (MYIMV tolerant) and VC1973A (MYIMV 
susceptible) (Ali et al. 1997), while NM10-12-1 is a mutant 
line from gamma irradiation whose original parent is not 
known. However, based on the pedigree provided by Ali et 
al. (1997), NM10-12-1 is possibly derived from 6601. 

In blackgram, resistance gene analog markers YR4 and 
CYR1 were reported to be completely linked with resistance 
to MYMIV (Maiti et al. 201 1). Marker CYR1 was also asso- 
ciated with resistance in mungbean (amplifies DNA band in 
resistant mungbean but not in susceptible mungbean). CYR1 
is proposed as part of the candidate disease resistance (R) 
gene (Maiti et al. 201 1). Recently, the R gene CYR1 was ful- 
ly isolated from blackgram (Maiti et al. 2012). Protein 
CYR1 may act as a signaling molecule to protect blackgram 
from MYMIV and be involved in recognizing the effector 
molecule of the pathosystem to contribute to incompatible 
interactions (Maiti et al. 2012). In the present study, markers 
YR4 and CYR1 were used in a polymorphism survey of the 
parents, but none showed polymorphism. Although CRY1 
showed no polymorphism between KPS2 and NM10-12-1, 
the association between CYR1 and MYMIV resistance in 
mungbean cannot be ruled out. The contrast between our re- 
sult and that of Maiti et al. (201 1) may possibly be due to the 
difference in genotypes of mungbean. DNA B of a MYMV 
isolate is an important determinant of host range between 
mungbean and blackgram. Balaji et al. (2004) showed that 
agroinoculation of mungbean and blackgram with two dif- 
ferent DNA B components together with DNA A of a black- 
gram isolate of MYMV- Vig resulted in distinctly different 
symptoms between the two crops. Additional research is re- 
quired to clarify whether CYR1 is involved in MYMIV resis- 
tance in mungbean. In addition, SSR marker CEDG180, 
which was reportedly associated with a major gene control- 
ling MYMIV resistance in blackgram (Gupta et al. 2013), 
showed no polymorphism between our mapping parents. 
Based on the linkage map of mungbean reported by Isemura 
et al. (2012), CEDG180 is mapped to LG10. This suggests 
that the genetics of resistance to MYMIV in mungbean and 
blackgram are different. 
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